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Abstract 

We address a conflicting report on the value and uncertainty of the astrophysical cross section factor of the *^C(a, 7 ) reaction 
extracted from existing data. In sharp contrast to previously reported ambiguities (by up to a factor 8 ), Schuermann et al. suggest 
an accuracy of 12%. We demonstrate that the so claimed “rigorous data selection criteria” used by Schuermann et al. relies on 
the s-factors extracted by Assuncao et al. But these results were shown in a later analysis (by this author) to have large error 
bars (considerably larger than claimed by Assuncao et al.) which render these data not appropriate for a rigorous analysis. When 
their “rigorous data selection” is adjusted to remove the results of Assuncao et al. the astrophysical cross section factor cannot 
be extracted with 12 % accuracy, or even close to it. Such data on the S e 2 values at low energies deviate by up to a factor two 
from their fit and exhibit a sharper slope rising toward low energies, leading to strong doubt on their extrapolated S £2(300) value 
and the quoted small error bar. Contrary to their claim the small value of S £i(300) ~ 10 keVb cannot be ruled out by current 
data including the most modem gamma-ray data. As previously observed by several authors current data reveal ambiguities in the 
value of S£i( 300) = approximately 10 keVb or approximately 80 keVb, and the new ambiguity that was recently revealed (by this 
author) of 5 £2 (300) = approximately 60 keVb or approximately 154 keVb, appear to be a more reasonable evaluation the status of 
current data. 
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During stellar helium burning oxygen is formed by the fusion of helium with carbon and is denoted by the *^C(a, 7 ) 
reaction. It has been announced three decades ago III] as the most irnportant and uncertain nuclear input of stellar 
evolution and it remains so today. In sharp contrast Schuermann et al. yfl claimed that the astrophysical cross section 
factor of the *^C(a, 7 ) reaction (as defined for example in |ll]) can be deduced with high accuracy of approximately 
12%, using a global R-matrix analysis of existing data. We demonstrate that existing data do not permit the claimed 
12% accuracy, or even close to it, not withstanding sophisticated R-matrix analyses. We also raise strong doubts on 
the values of the astrophysical cross section factors quoted by Schuermann et al. 12] • 

Two comments are in order from the outset. First, the “rigorous data selection criteria” of Schuermann et al. 12] 
include ten data points shown with the measured 100% (or more) error bars. Such data points should be considered 
as upper limits and as such they cannot be included in a rigorous chi-square analysis, since the contribution of an 
upper limit to chi-square cannot be rigorously evaluated. Most bothersome are the five such data points shown at 
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low energies below 2 MeV (2 data points for S e\ and 3 for Sei)- More significant as we discuss below, other data 
included by Schuermann et al. ||2l were shown in a later analysis (by this author fl) to have similar large error bars 
(close to 100%) which also render these data not useful for a rigorous chi-sqaure analysis. 

Second, the data points below 2 MeV are essential for having high sensitivity to the contributions of the bound 1“ 
and 2^ states at 7.1169 and 6.9171MeV in '®0, respectively, that govern the value of the astrophysical cross section 
factors at stellar energies: 5^1 (300) and 5^2(300), respectively. But the contribution to^^ of these low energy data 
points is overwhelmed by the large number of high energy data points included by Schuermann et al, leading to a 
reduced sensitivity to the low energy data. In addition the inclusion of higher energy data points raises additional 
question(s): for example concerning the sign of the interference of the 2^ at Ecm = 2.68 MeV that was shown to be 
positive and leading to an extrapolated 5^2(300) = 62^g keVb fllH. This S-factor is 16% lower than quoted by 
Schuermann et al. with a difference that is considerably larger than the quoted uncertainty of ^£2(300). We note 
that the exclusion of high energy data points and consideration of data points only below 1.7 MeV was shown by this 
author to lead to very different conclusions as we discuss below. 

More troubling is the inclusion and in fact their reliance on the published data of Assuncao et al. |@1- We note 
from the outset that Brune et al. l|3l concluded several systematical problem with the data of Assuncao et al. |6]. 
These data |0 were re-analyzed by this author (|3l and I revealed very large error bars of the extracted s-factors |3], 
considerably larger than stated by Assuncao et al. |@]. 

We recapitulate here a few observation made in my chi-square analysis J^l of the data of Assuncao et al. 1@]. For 
example even though no discernible peaks can be established in the gamma-ray spectra measured at 90° - 130° shown 
in Fig. 6 of Assuncao et al. |@], an angular distribution is claimed to have been measured at these backward angles 
at the indicated energy of Eajat - 1 -850 MeV. I showed that the angular distribution labeled as Ea,iab - 1 -900 MeV 
(E = 1.340 MeV) |@1 can be fitted with values that vary by a factor of six with similar reducedvalues; see Fig. 
1 of Ref. 1^. My chi-square analysis of all data points measured by Assuncao et al. |@1 below 1.7 MeV lead to 
values that vary by a large factor without significant variation in chi-square. The resultanting uncertainties are shown 
in Fig. 2 of my paper J^l and they are close to 100%. As discussed above such data points must be considered as upper 
limits and the data of Assuncao et al. |@] cannot be included in a rigorous chi-square analysis. We emphasize that 
excluding the results of Assuncao et al. is not a matter of choice for “data selection”, rather it is dictate by standard 
considerations of a rigorous chi-square analysis. 

In the same paper 1 also demonstrated that the disagreement of the E1-E2 measured phase angle (^ 12 ) with (the¬ 
oretical) predictions, shown in Eig. 11 of Assuncao et al. ||3l, is a violation of unitarity ||3l and not just a mere 
disagreement with the prediction of R-matrix theory as suggested by them |@]. A violation of unitarity is a clear 
indication of serious problems with the data or the data analysis. 

We conclude that the “rigorous data selection” employed by Schuermann et al. (|2l should not include the data of 
Assuncao et al. and we remove the data of Assuncao et al. from the data satmle analyzed by Schuermann et al. 
But since they already removed the data of Redder et al. and Ouellet et al. H] from their choice for the current 
precise data, the adjusted choice of data of Sei measured below 2 MeV includes only the data of Kunz et al. ifTol] . 
This is clearly an unacceptable situation that must be alleviated by new data measured at low energies and it cannot 
be remedied by sophisticated global R-matrix analyses. In the same time it is important to comment here that my 
analysis of the data of of Kunz et al. Jl^l agrees with Kunz et al. as stated in my paper iQ]. 

Nevertheless the comparison of the low energy R-matrix curve shown in Eig. 1, taken from their Eig. 4 ||2l, reveals 
a disagreement by up to a factor of 2 with the adjusted choice of data sample used by Schuermann et al. (which does 
not include the data of Assuncao et al. and only includes the data of Kunz et al., excluding the data points measured 
with 100% error bars). More troubling is the fact that in comparison to the data their R-matrix fit has the wrong slope 
at low energies as shown in Eig. 1. Clearly these data points below 2 MeV are crucial for delineating the contribution 
of the bound 2^ state of and the slope is very important for an accurate extrapolation of S £2(300). 

We conclude that when considering the adjusted choice of the data sample used by Schuermann et al., the value 
they quote for S e 2 cannot be substantiated and certainly we cannot support their claimed small uncertainty of the 
extracted ^£2(300). We refer the reader to Ref. I^l for a complete chi-square analysis of all currently available data 
below 1.7 MeV that bifurcate and lead to ambiguity in the extracted value of S £2(300) = 60 + 12 or 154 + 31 keVb. 

Concerning the value of 5£i(300): Schuermann et al. (|2l consider the possibility of destructive interference of 
the K sub-threshold state at 7.1169 MeV and the 2.42 MeV K resonance in and they claim “the constructive 
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Figure 1. (Color online) Comparison of the low energy portion of the R-matrix fit of Schuermann et al. J3] with the S e 2 data measured by Kunz 
et al. ca. As discussed in the text the data of Assuncao et al. H have been removed from their “rigorously selected data” leading to a major 
discrepancy with the remainder of the data (i.e. the Kunz et al. data) and specifically with the wrong predicted slope. 

solution is strongly favored and the destructive interference pattern has been rejected”. This claim is based on the 
obtained= 265 and 233 for 5£i(300) = 7.9 and 83.4 keVb, respectively. When considering the large number of 
capture data points included by Schuermann et al. (243) we conclude that the obtained difference of 32 is hardly 
significant to warrant the rejection of the destructive interference pattern [yielding the small 5£i (300) = 7.9 keVb]. 

Furthermore, since they state that a one sigma variation of the total s-factor leads to Ax^ 21 [121], the obtained 
X^ difference of 32 is not significant enough and certainly far from the usual five sigma used to substantiate a claim. 
Clearly both fits have reduced values that are close to unity, and the fit that leads to the small 5^1 (300) = 7.9 
keVb cannot be considered as ruled out by the data since the reduced x^ is close to unity. As such the destructive 
interference cannot be ruled out with the certainty claimed by Schuermann et al. It is important to note that it is not 
sufficient to demonstrate the good fit for the large S ei (300) solution, but one must also rule out the small S ei (300) 
solution. This has not been achieved by Schuermann et al. and current data do not allow us to rule out the small 
5^1(300) solution, leading to the ambiguous value of 5^1(300) [jit]. 

The rejection of the small S ei solution by Schuermann et al. based on x^ consideration is made more doubtful 
since they included and relied on the data of Assuncao et al. which were found to have unrealistic small error bars 
This makes the minuscule difference in the reduced;i^^ even more troublesome for rejecting the destructive interference 
pattern. 

Similarly the fit to the modern data that was published in a peer reviewed conference proceedings iH (and was 
neglected by Schuermann et al.) states the numerical values: 5^1 (300) = 77.9 and 4.3 keVb with = 9.0 and 9.6, 
respectively, see Fig. 5 of lEl. Such a small difference in x^ in of itself does not allow rejecting the small value 
solution and the need to re-evaluate the error bars quoted in Ref. |3, 11| weakens the possibility of rejecting the small 
s-factor solution using these modern data alone. We thus conclude that this modern gamma-ray data analysis id just 


the same as the analysis of Schuermann et al. (|2l and previous data analyses 112, 13] do not support ruling out the 
small 5 £1 ss lOkeVb solution. 

We conclude that a realistic evaluation of current data does not permit the determination of the astrophysical cross 
section factors with the 12% accuracy suggested by Schuermann et al., nor can we support their claimed values of 
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5£i( 300) and SeiOOO)- A more suitable conclusion is that both extrapolated SeiQOQ) and 5£2(300) cross section 
factors are ambiguous with the values listed in my previous publication |13l- 

Unlike the strong claim of 12% accuracy suggested by Schuermann et al. the observed large ambiguities justify 
and promote a new and different research effort to determine the astrophysical cross section factors with the required 
uncertainty of 10% or better. Indeed proposals for determining the cross section at very low energy have been de¬ 
veloped for the HlyS gamma-ray facility in the USA Q, and the newly constructed ELI-NP facility in Bucharest 
as shown in Q. Measurements with gamma-beams are favored by the detailed balance factors of 50 - 100 and are 
made possible due to the anticipated high intensity (lO^y/i) that will allow a measurement at E^.m. = 1.0 MeV within 
two weeks of beam time. The design goal sensitivity of these measurements which is shown in Eigs. 2 and 3 of lll5n . 
promises to resolve the observed ambiguities in 5£i(300) and 5 £2 (300). 

In closing we note that for example as shown in Eig. 1 of jlbll the suggested value of the total astrophysical cross 
section factor quoted in ^ of 161 + 19-1-8-2 keVb (i.e. a multiplicative factor of 1.61 as defined in 116]) is exactly 
at the boundary (170 keVb) where a 25 solar masses star is predicted to be oxygen rich (C/O < 1) and thus skip the 
carbon burning stage and collapse to a black hole. Thus a resolution of the ambiguities in S £i (300) (approximately 
10 or 80 keVb) and ^£2(300) (approximately 60 or 154 keVb) noted in Ref 13] is essential for progress in stellar 
evolution theory. 
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